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Abstract 

The potentials between two .B-mesons are computed in the heavy-quark limit using quenched lattice 
QCD at m n ~ 400 MeV. Non-zero central potentials are clearly evident in all four spin-isospin 
channels, (I, si) = (0, 0), (0, 1), (1, 0), (1, 1), where s/ is the total spin of the light degrees of freedom. 
At short distance, we find repulsion in the I ^ si channels and attraction in the I = si channels. 
Linear combinations of these potentials that have well-defined spin and isospin in the i-channel are 
found, in three of the four cases, to have substantially smaller uncertainties than the potentials 
defined with the s-channel (I,S[), and allow quenching artifacts from single hairpin exchange to 
be isolated. The BB*ir coupling extracted from the long-distance behavior of the finite-volume 
t-channel potential is found to be consistent with quenched calculations of the matrix element of 
the isovector axial-current. The tensor potentials in both of the s/ = 1 channels are found to be 
consistent with zero within calculational uncertainties. 
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I. INTRODUCTION 



Nuclei and nuclear processes can be described with remarkable precision by treating the 
nucleons as non-relativistic particles interacting via a local potential. The wealth of nucleon- 
nucleon (NN) scattering data has enabled the construction of precise phenomenological 
potentials, defined up to unitary transformations, which are used to calculate the spectra 
of the light nuclei. As expected from QCD, the two-body potentials alone are insufficient 
to reproduce the spectra, but when supplemented with three- and four- (and higher) body 
interactions successfully reproduce the structure of light nuclei [l|, 0] • During the last 15 years 
or so, immense effort has been put into developing an effective field theory (EFT) to describe 
the interactions between nucleons [3,4], and to allow for a systematic improvement of nuclear 
physics phenomenology using the approximate chiral symmetry of QCD. The values of the 
counterterms that appear at a given order in the EFT expansion have to be obtained from 
experiment or lattice QCD. In situations where experiments are not possible, lattice QCD 
is the only rigorous calculational technique with which to determine such counterterms. 
Recently, NN scattering (EJ and hyperon-nucleon scattering 0] have been calculated with 
fully-dynamical lattice QCD 1 by measuring the finite-volume shifts of two particle energies 
0, H[ • Due to limited computational resources, the calculations were performed at unphysical 
pion masses, 350 MeV, at the upper limits of the range of applicability of the EFTs. 

Until the computational resources for lattice QCD calculations are significantly greater than 
presently available, it will not be possible to calculate NN scattering parameters at a large 
number of different energies and then construct NN potentials in the same way that NN cross- 
section measurements are processed. However, in addition to extracting the NN phase-shifts 
at various moment, one can hope to learn qualitative information about the EFT describing 
NN interactions by performing lattice QCD calculations of systems that are similar. 

In this work we study the potential between two mesons in the heavy-quark limit [1, 
10, 11, 13] 2 , a limit in which the potential is a well-defined object. .B-mesons are isospin- 



| hadrons, and in the heavy-quark limit the spin of the light degrees of freedom (ldof) 
becomes a good quantum number, si — ~, as spin-dependent interactions with the heavy- 
quark are suppressed by l/mi,. At distances which are large compared with the chiral 
symmetry breaking scale, A x , the EFT describing the interactions between two .B-mesons 
is the same as that between two nucleons as the isospin-spin quantum numbers are the 
same. The differences between the two EFTs are in the values of the counterterms. At 
distances that are of order, or shorter than A x the interactions between two .B-mesons will 
be arbitrarily different from that between two-nucleons as the structure of the hadrons are 
very different, and in particular that strong-Coulomb interaction between the heavy-quarks 
becomes dominant, behaving as ~ a(r _1 )/r as the separation, r — > 0. 

In addition to providing insight into the NN potential, the potentials between B-mesons 
are interesting in their own right. A precise determination of the potentials between two 
B-mesons will allow for investigations of possible shallow bound states. These would be 
molecular tetra-quark states, similar to the deuteron in the NN-sector. The location of such 



1 Recently, it has been claimed that the nucleon-nucleon potential has been calculated with quenched lattice 
QCD. We believe that the arguments used to extract the potential are flawed, and discuss this further in 
Appendix [A~J 

2 The B and B* mesons are degenerate in the heavy-quark limit and henceforth we will use B-meson to 
denote the B and B* super-multiplet. 



2 



molecular states would be very sensitive to the potentials (due to the fine-tunings) and as 
such, quenched calculations at the unphysical pion mass would in general provide unreliable 
results. 

Lattice calculations of the potentials between two S-mesons in the heavy-quark limit 



have been performed previously [13j, [14|, [15j, [16|, [17|, UM, uM, 120, [21|, [22], |23J]. However, 



given the large statistical uncertainties in those calculations, the potentials remain largely 
unexplored. We have chosen to work in a relatively small lattice volume in order to explore 
the intermediate and short- distance components of the potential, but this is at the expense 
of having contributions from image pairs somewhat mask the long-distance component. 
Further, we have attempted to extract the tensor potentials in the si = 1 channels, but have 
not found results that are statistically different from zero. Due to limited computational 
resources, we have performed quenched calculations. It is important to stress that the long- 
distance component of the potential computed in quenched QCD is polluted by the presence 
of "one-hairpin-exchange" (OHE), as discussed in Ref. (24J, which becomes dominant at large 
distances due to its exponential fall-off, as opposed to the Yukawa-type behavior of one pion 
exchange (OPE). However, the OHE contribution can be isolated by defining potentials with 
well-defined t-channel spin-isospin quantum numbers. In three of the t-channel potentials, 
quenching artifacts are expected to be higher order in the quenched EFT expansion. Using 
these finite-volume t-channel potentials, we are able to investigate the long range part of 
the infinite-volume potential and extract a BB*n coupling consistent with that measured in 
quenched lattice calculations of axial matrix elements. 

The outline of our paper is as follows. In Section [III we discuss the numerical implemen- 
tation of our calculation. Sections II I II and HVl present our single particle results for the heavy 
hadron spectrum, including exotic states. Section |V] presents the main results of our work, 
the BB potentials in the various channels. These results are then discussed in Section IVII 
In Appendix [Aj we discuss the issue of extraction of nucleon potentials from lattice QCD 
wavef unctions, and in Appendix [B] we provide details of the perturbative lattice calculations 
needed in this work. 



II. DETAILS OF THE LATTICE CALCULATION 

Our calculations were performed using 284 quenched configurations of dimension 16 3 x 32 
generated with the DBW2 action [25|, 26[ for (3 = 1.04, giving a lattice spacing of 6 = 



0.0997 ± 0.0015 [27]. On each of these configurations, eight Wilson light-quark propagators, 
equally spaced in the time-direction and offset in space, were generated from smeared sources 
to determine the light hadron spectrum. The light-quark mass selected gave rise to a pion 
mass of m n = 402.5 ±6.7 MeV, and the other hadron masses that are shown in Table. [B 
The finite lattice spacing and finite-volume effects [28| have not been removed from these 
masses and are expected to be a few percent. 

In order to compute the potential between two fi-mesons separated by lattice vectors 
r = b(n, 0, 0), 6(0, 0, n) for n = 0, . . . , 8 and additionally by r = 6(1, 1, 0), 6(2, 1, 0), Wilson 
light-quark propagators were generated from smeared sources on one time-slice located at 
each point in the x-z plane on two adjacent spatial-slices in the y-direction on each gauge 
configuration. Therefore, a total of ~ 1.4 x 10 5 light-quark propagators were generated. 
This choice of lattice separation vectors was dictated by the available computational time, 
and not by physics. The calculations were performed on a 16-node dual-Xeon cluster and a 
number of workstations. The total computational cost of this work was ~ 40 Gflop-yrs. 
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Quantity b M M [MeV] 



m N 
rriA 



0.2034 ± 0.0015 
0.3754 ± 0.0080 
0.5756 ± 0.0080 
0.6770 ± 0.0095 
0.102 ±0.014 



402.5 ±6.7 
743 ± 19 
1139 ± 23 
1340 ± 28 
201 ± 27 



TABLE I: The masses of the light hadrons at finite lattice spacing and finite volume. 

In the heavy-quark limit, the heavy-quark propagator is the tensor product of a Wilson 
line and a positive-energy projector, 



S Q (x,t;t ) = 



74 



n tMx,f), 



(i) 



t'=t 



where U^(x) are the gauge link variables and the product is time-ordered. Our Dirac matrices 
use the Euclidean Dirac convention. The light quark propagator, S(x.,t;y,to), is generated 
with the unimproved Wilson action, thereby introducing 0{b) discretisation errors. It is 
generated from a gauge-invariant Gaussian smeared-source. 

To determine the single particle energies of the heavy hadrons, the correlators 

Cb(Mo) = E tr [^(x,t;to)5 t (x,t;x,to)] #(x) , (2) 

x 

C S6/Ab (t, to) = E C( x ' *; to)e ijk e^' k ' (s«'(x, t; x, t )v) (rS"'(x, t; x, to)) H(x) , 



(3) 



were computed, where the Dirac matrices are T = C75, C'-fi for the and S| respectively. 
The trace is over color and spinor indices, and the function H(x) is unity if a heavy-quark 
source was placed at the point x, and vanishes elsewhere. In the baryon correlators, upper 
indices label color and lower Greek indices label spin. 

In order to measure the potential, we computed the correlators Ci tSl (t,to,r) given by 



C ,o(Mo,r) 
Ci )0 (Mo,f) 

Cb,i(Mo,f) 

Ci,i(*,t ,f) 



£ K + A 



ff(x), 



+ A a) + A 



(2) 



A (2)' 



#(x) 
HU) 



(4) 
(5) 

(6) 

(7) 



where 



A { o ] = tr D [tr c (5 t (x,t;x±r,t ) 1 SQ(x,t;t ))r 

xtr c [S {^(x + r, t; x, t )S Q (x + r, i; *„)} 



(8) 



(a) 



(b) 



FIG. 1: Disconnected and connected quark contractions contributing to the correlation functions 
in Eqs. ([H)-©. Heavy and light lines correspond to the heavy- and light-quark propagators 
respectively. 
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(2) 



-tr 



S |S ,t (x,t;x + r,t )5'Q(x + r,t;to)} 



A 



(i) 

± 



tr 



xr 

r ± Sg(x,t;t )S ,t (x,t;x,to 
xtr 



S t (x + r,t;x,t )l S Q (x,t;to)T. 



(9) 
(10) 



1 ± 



-tr 



T±S {S f (x + r, t; x + r, t )S Q (x + r, t; to)}] , 
Sq(x, t; t )S {^(x, t; x + r, t Q )S Q (x + r, i; *„)} T±S*(x + r, i; x, t )T± , (11) 



T = ^(1 + 74)75 and T± = |(1 ±z'7i7 2 )(l +74) and tr c , tr D and tr indicate traces over color, 
spin and both. Here S{. . .} indicates translation of the propagators by the lattice vector — r 
and the transpose T denotes spin transpose only. The different contributions to the various 
correlators, A^\ and A^\_ correspond to the two contractions shown in Fig. [U 

For each correlator, we determine the ground state energy by seeking plateaus in the 
ensemble jackknife average of the effective energy, 



6 *o) = log 



Ci, S[ (t — I? to) 

Cl,si(t, to) 



(12) 



III. THE HEAVY HADRON SPECTRUM 



In the heavy quark limit, the mass of the B meson is 



M. 



B 



m b 



An + 0(l/m b ) 



(13) 



where m b is the heavy quark mass and Aj tSl denotes the energy of the ldof with total isospin 

/ and spin si. To determine Ai 1 from lattice calculations, the energy, £11, of a meson 

2 ' 2 2 ' 2 

composed of a Wilson-line (static color source) and a light anti-quark is computed using 
Eqs. (j2j) and f[T2l . This by itself does not isolate An, as the interactions of the static 



source with the gauge fields generate a residual mass |29j for the heavy-quark, 5m, which 
while vanishing in dimensional regularization, is non-zero on the lattice and scales as 1/6 
with the lattice spacing 3(|. Therefore, both Si 1 and 5m diverge as 1/b, but the difference 



between them is finite in the continuum limit, and is Ai 1 . 

2'2 



Ai 1 



Si 1 

2 ' 2 



5m 



(14) 
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and more generally, A/ iSj = Si tSl — 8m. 

The residual mass of the static source has been computed previously out to the two-loop 
level in quenched lattice perturbation theory for the Wilson action [31]. At the one loop 
level, the residual mass for a gauge action, /, is given by 



dq x I dq y I dq z G$ (q x , %, ?*, 0) 



qb is dimensionless. G^(k 



(15) 



where g 7 = 2 sin (<? 7 /2), and q = qb is dimensionless. G)^{k x , k y , k z , k t ) is the lattice gluon 
propagator for the particular gauge action, /. This has the form G^ llsm \k x ,k y ,k z ,kt) 
Ijk 2 for the Wilson action, but is considerably more complicated for improved actions, such 
as the Liischer-Weisz (LW) [ 



33 



26 , and Iwasaki \3m actions. For these 



DBW2 [25|, 

actions, the form of the propagator was presented in Ref. [35| and involves an improvement 
coefficient, c u with c LW = -1/12, cf BW2 = -1.40686, and 4 wasaki = -0.331 (for a recent 
review of lattice perturbation theory, see Ref. [36j]). The quantity has been computed 
previously for the Wilson action, ^C 1 ^ 4 ^ ™) = 2.1173, and for the improved actions we find 
that S^ DBW2 ^ = 0.6921, S^ LW ^ = 1.8335 and S^ Iwasaki ^ = 1.3598 (these numbers differ from 
those presented in Ref. [37| where an incorrect improved gluon propagator was used). At 
the one-loop-level, the choice of scale is not well-defined, but as the only scale in the lattice 
calculation is the lattice spacing, it is convenient to use /i = 1/b. 3 Therefore, at the one-loop 
level, and using 5(2 GeV) = 0.299 ± 0.015 4 , we have 5m££ W2 = 410 ± 21 MeV for the 
DBW2 action. 

One can make an improved estimate of the residual mass term by using the Brodsky- 



Lepage-Mackenzie (BLM) scale-setting procedure [38(|, which includes the part of the two 



loop contribution arising from the running of the strong coupling over the momenta in the 
one-loop diagram. This BLM improved residual mass is given by 5 



S (f) 



a 2 (/j,)f3 
Anb 



T (/) (/i), 



where 



T u \v) = I d ^ x d Qy / d & G oo 0) log 

07T J —TV J —TV J —TV 



f x + € + fz 



fi 2 b 2 



(16) 



(17) 



and (3 — 11 for quenched QCD. The a 2 f3 term in Eq. (fl~6|) can be perturbatively removed by 
defining the BLM scale, <f (/) , such that log (V (/) &) = T (/) (6 _1 )/(25 (/) ), and therefore, 



S if) 



(18) 



3 Perhaps a better estimate would be fi = ir/b as that is the maximum momentum in the one- loop diagram. 
The typical momenta in the one- loop diagram to be somewhat less than this. 

4 The strong coupling on these DBW2 lattices has been determined to be = 0.154 in the MS- 
scheme j^l, significantly smaller than the experimentally constrained value of a(2 GeV) = 0.299 ±0.015. 
This suggest that the perturbative relation between the value of the plaquette and the strong coupling is 
only slowly convergent. 

5 Our definition of the improved residual mass is different than that in Ref. [30( but agrees with Ref. [31 1 . 
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FIG. 2: Effective mass ratios of the correlation functions for the B (red), A& (green) and £& 
(violet). The shaded regions correspond to the fit range, central value and uncertainty (statistical 
and systematic errors are added in quadrature) for each hadron. 



For the Wilson action, we find 'X( Wllson )[p x ) — 1.562 which produces a BLM scale 
q*(Wiison) _ i 445/6 that has previously been shown to accurately estimate the full two- 



loop result for the residual mass using the MS coupling [3l|. For the improved actions, we 



find that r( DW2 )(6" 1 ) = -0.239, which leads to q < DBW2 ) = 0.841/6, T^ LW \b^) = 1.121, 
which leads to q* {LW) = 1.358/6 and for the Iwasaki action, T ( Iwasaki ) (fe- 1 ) = 0.437, which 
leads to q*V™ asaki ) = 1.174/6. Therefore, using a(q< DBW2 *>) = 0.326 ± 0.018, the BLM 

improved estimate of the residual mass for our DBW2 lattices is Sm^T = 447 ±25 MeV. 

In Table [TT] we present the extracted lattice energies and resultant energies of the ldof 
for the .B-meson and the and E& heavy baryons. The effective mass ratios corresponding 
to these measurements are shown in Fig. [2j We have performed both correlated and un- 



6 At finite lattice spacing, the definition of the BLM scale becomes ambiguous as the two loop contributions 
that the BLM procedure is attempting to resum become dependent on the details of the discretisation. 
In particular the continuum log(|<f| 2 ) becomes a complicated function of the lattice momenta and im- 
provement coefficients. A full two-loop calculation will be required to determine the efficacy of our BLM 
estimate using a(q*) for the improved actions and a priori there is no reason to assume the agreement 
found for the unimproved Wilson action persists in these cases. As an indication of possible lattice ar- 
tifacts in the definition of T we have also computed q*^ by replacing log(|g| 2 ) — > log(|<7| 2 ) in Eq. [fTTj). 
finding q<wiison) = 1.671/6, q<DBW2) = 0.918/&, = 1.552/6, and q< Iwasaki ) = 1.318/6. Similar 

perturbative shifts in the BLM scale are induced by self-consistently evaluating T at the scale q*v>. 



7 



correlated single and double exponential fits to the correlation functions in order to extract 
the ground state energy and to estimate the uncertainty in the extraction. The differences 
between the extracted ground state energies from these procedures is encapsulated in the 
systematic error. The fitting range quoted in Table [II] (and all other tables) is that used in 
fitting a single exponential to the correlation function. The statistical errors are determined 
by the Jackknife procedure, omitting a single configuration at each evaluation. Further, a 
Bootstrap analysis was also performed on the data, with both techniques providing similar 
central values and uncertainties. To eliminate uncertainties common to the energies of two 
different hadrons, we formed the correlated differences between the energies of the ldof in 
the various systems. The residual masses of the static sources cancel in these combinations 
and the results are displayed in Table IIII1 



Hadron fit range b £(DBW2) £ (DBW2) [MeV] A^bwP [MeV] 



B 


11 -> 15 


0.5539(37)(50) 


1096(18)(10) 


649(31)(10) 




9 -» 14 


0.7934(65)(22) 


1570(27) (04) 


1123(36) (04) 


S 6 


10 -> 14 


0.8575(74) (74) 


1697(29)(15) 


1250(38)(15) 



TABLE II: The spectrum of hadrons comprised of one static color 3 source and light (anti-)quarks, 
for a light quark mass giving the light hadron spectrum shown in Table. [H g( DBW2 ) [$ the mass 
of the hadron determined in the lattice calculation. ^dbW'P = 8 <yDBW2 ^ — Sm^g^-^ ^ s the BLM- 
improved determination of the energy of the ldof. The first uncertainty is statistical and the second 
uncertainty is systematic. 



Hadrons bS = bSA SA [MeV] 



A b -B 


0.2395(75)(72) 


474(15)(14) 


Z b -B 


0.304(08)(12) 


601(16)(24) 


Z b -A b 


0.0641(98)(96) 


126(19)(19) 



TABLE III: The mass differences between the hadrons comprised of one static color 3 source and 
light (anti-)quarks. The residual mass of the static source cancels in the differences, leaving the 
differences between the energy of the light degrees of freedom, denoted by 5A. The first uncertainty 
is statistical and the second is systematic. 



IV. EXOTIC BARYONS 



As a by-product of computing the potential between two 5-mesons, we have computed 
the masses of baryons formed from a static-source transforming in the 6 of color and two 
light anti-quarks with (I,si) = (0,1), which we denote as Ag, or (J, Sj) = (1,0), which we 
denote as Eg. These are calculated by putting two static sources, each in the 3 of color at 
the same point in space, and requiring the light anti-quarks to have the appropriate values 
of (I, si). As the spin of the ldof decouples from the static source(s) and is a good quantum 
number, the ldof form a color 6. While the Casimir of the 3 representation is C(3) = 4/3, 



S 



1.4 



1.3 



bE 1.2 



* 

^6 



1.1 



10 



12 



t/b 



FIG. 3: The effective mass plots for the exotic hadrons, Ag and Eg. 



Hadron 


flt 


range 


6 g(DBW2) 


£ (DBW2) [ Me y 


, x («+« 2 /3) [M V1 

J 6;DBPV2 I ev J 




6 


-> 10 


1.253(08)(02) 


2480(16) (04) 


1364(64) (04)) 


^6 


7 


-> 11 


1.278(10)(05) 


2529(20)(10) 


1413(65)(10) 












5A 








0.025(13)(05) 




49(25)(11) 



TABLE IV: The spectrum of exotic baryons comprised of a static color 6 source and light anti- 
quarks, for a light quark mass giving the light hadron spectrum shown in Table. HI £( DBW2 ) [ s 
the mass of the hadron determined in the lattice calculation. Using the BLM- improved one-loop 
lattice perturbation theory calculation of the residual mass, as discussed in the text, we present 

where the residual mass term is evaluated at the 



results for A^ a Jl = £^ DBW2 ^ - 5m {a+a2p) 



— (a+a 2 f: 

^6;DBW2 ~ " "■' l 6;DBW2' 

scale q*( DBW2 ) given in the text. The difference between the masses of the exotic states is given 
in the last line. 



it is C(6) = 10/3 for the 6 representation, and the residual mass that must be subtracted 

from the energy calculated on the lattice is Sm^swl = 1116 ± 62 MeV. The effective 
mass plots for these hadrons are shown in Fig. [3] and the energies of the ldof are shown 
in Table IIVI These states are considerably more massive than the non-exotic hadrons and 
achieve plateaus at earlier times. As such exotic baryons have not been observed and heavy 
quarks transforming in the 6 (or other higher representations) of color are not required in 
nature (see Refs. jH, 0, BlU for a number of proposals) and are not expected to be found, 
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we do not dwell on these results further. 



V. B-MESON POTENTIALS 



In the heavy-quark limit, the separation between the two B-mesons is a good quantum 
number and the potential, V(r), is simply defined by the difference between the energy of 
the two S-mesons separated by a displacement vector r (defined by the separation of the 
static 3 color sources) and the energy of two infinitely separated 5-mesons, 

V I<St (r) = E I)Sl (v) - 2M B = A M (r) - 2A |; i , (19) 

where the isospin and spin of the ldof can take the values I = 0, 1 and si = 0, l. 7 Contributing 
to both of these energies are the interactions between the ldof, the interactions between the 
ldof and the static sources, and the interactions between the static sources. As |r| — > oo, each 
of the contributions factorize, leaving the contribution from two non-interacting 5-mesons, 
and therefore a vanishing potential. 

The lattice calculation is a straightforward extension of the calculation of A for the heavy 
hadrons. The energy, £ff(r), of two B-mesons composed of static sources and light-quarks, 
placed on the lattice with relative displacement r, is computed using the correlators in 
Section [Til The potential computed on the lattice then becomes 

VS?(r) = £%{r) - 2 £ U , (20) 

where 



£ff l (v)=A l!Sl (v)-8V R (r) + 28m , (21) 

(the subscript R labels the color representation of the heavy quark system dictated by the 
light quantum numbers I and si). The residual masses of the static sources induced by 
interactions with the gauge fields cancel in Vj^*(r). However, a perturbative subtraction 
corresponding to the differences in interactions between the static sources in the continuum 
and on the lattice, (which at leading order in the strong coupling arises from one gluon 
exchange (OGE)) 8Vr(t) = Vqge — "^oge remains. Therefore, the energies measured in the 
lattice calculation, and the potential between two 5-mesons in the continuum are related 
by 

Vr f-l (r) = £f£(r) - 2£^ + 8V R (r) . (22) 

In the cases where the spin of the ldof is si = 1, the potential receives contributions from 
both a central component and a tensor component, 

V^ s=1 \r) = Vc(r) + S 12 V T (r) , (23) 



7 We classify these states using the quantum numbers of the infinite-volume continuum. For s; = 0, 1 this 
is legitimate as there is a exact correspondence between these representations of 0(3) and the A\ and 7\ 
representations of H(3) 42[. The continuum symmetry is 0(3) as the presence of the infinitely massive 
quarks breaks 0(4) to its spatial subgroup. 



10 



where 

5i2 = ^(S + (f x ~tr y ) + S4f x + tf y ) + 2S z f z ) 2 - 2S 2 , (24) 

Si are the spin operators, and r is the unit vector in the direction of the displacement. It 
follows that the central and tensor potentials can be determined from the potential at x and 
z displacements, 

Vt l \r) = ~ (> (S=1 V<y + 2^)^,)) 

vf =1) (r) = \ ( V^\r e z ) - V^\r e x ) ) , (25) 
where ej is the unit-vector in the "j" direction. 

A. Lattice Spacing Effects in a Finite- Volume 

The potential measured in the lattice calculation will differ from that at infinite-volume 
due to the presence of image 5-mesons resulting from the periodic boundary conditions in 
the spatial directions of the lattice. Therefore, the single particle energies that are extracted 
correspond to the energy of single particle that is interacting with its images, located at 
|r| > L. In the case of the energy of two particles interacting in a periodic cubic volume, 
the potential energy, V^ L \r), measured includes the sum over the contributions from the 
images 8 

V (L) (r) = V(r) + V{y + iiL). (26) 

n^O 

When the displacement between the mesons is |r| > L/2 the interaction with the nearest 
image is more important than the interaction within the volume. Consequently, we have 
only computed the potential for |r| < 8 lattice spacings on 16 3 x 32 lattices. After taking 
the continuum limit of the lattice calculation, the finite-volume effects due to the images 
must be removed to recover the infinite- volume continuum limit potential, V(r). This is 
discussed below. 

The finite lattice spacing, b, eliminates ultraviolet modes on the lattice leaving |q| < n/b, 
and hence the strong Coulomb potential that exists between two static- sources due to OGE 
is significantly modified for |r| less than a few b. In the heavy-quark limit, the OGE potential 
is spin-independent but does depend upon the color representation of the combined heavy 
quark system. 

The potential between two static color 3 sources combined into a color 3, at a finite 
lattice spacing, b, and in a finite-volume becomes 

f dq x f dq y f c%G 00 (q)e^ f 

37TO J-1T J-TT J-n 

(f ) 3 "?* e ' 2m '" i Gw ( 2 sil1 (tO • 2 sin (x) ' 2 sin (x ) ' °) (27) 

This is correct for interactions via single particle exchange but receives corrections that we discuss in 
Section |VD] 



ylatt / \ 



«(/i) 

3^6 
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where L = L/b is the spatial extent of the lattice in lattice units, and f = r/b is the 
displacement between the static sources in lattice units. The summation in Eq. (J271) is over 
all -L/2 < n : 



xi lb yi 



< L/2. This finite- volume expression has an infrared divergence due 
to the n = mode, however, the difference between the finite-volume OGE potentials at 
finite lattice spacing and in the continuum is of the form 



5V^ L) (r) 



/2tt n 3 
3tt 2 6 VT 




E 



gi2-7rn-r/L 



Kl<4 



E i2Trn-r/L s-t 



oo (n x ,hy,h z ,0) 



(28) 



where hi = 2 sin (jnii/Lj, and is well-behaved in the infrared. Spurious contributions from 
ill-defined low-momentum gluon modes included in the above sums cancel to a large extent 
with residual effects at most of 0(b 2 ). Further discussion of this issue and the numerical 
evaluation of 5V^ lq {t) can be found in Appendix IB1 

The BLM procedure is again used to set the scale of the correction factor, leading to 



5Vt L) (r) 



a(/i) f2n 



' L\ ~ e a*n*/L f aU)B /47r 2 (|n| : 



3vr 2 6 V L 

- e^/ L G 00 (h x ,h y ,h z ,0)(l - ^^\og nil 



47T 



fl 2 b 2 



Arc 



a(q*(T,L)) 



A(t,L) + ... . 



(29) 



The coefficient functions A(r, L) and B(r , L, /x) evaluated on DBW2 lattices, using the tech- 
niques in Appendix [Bj are given in Table. |V] at the required separations. The resulting 
BLM-scale and potential shifts are given in Table. IVII (results for the color 6 OGE potential 



are related by SV^ N ^ LO (r) = — 1/2 5V^^ LO (r)). The correction factors, SV^ Q (r) and 
L ■ 1 should be added to the lattice measurements of Si s , — 26 1 i to give the potential. 

2 ' 2 

At relative displacements that are large compared with the lattice spacing, this correction 
factor scales as 6 2 /|r| 3 , as expected. However, there are still 0(b) lattice artifacts from the 
discretisation of the light quark sector. These can only be eliminated using data at different 
lattice spacings or using a light-quark action that is (9(6)-improved. 



,(L) 



(L) 



B. The Lattice and Continuum Finite- Volume Potentials 

Using the techniques described in Sec. HI] we have computed the correlation functions 
corresponding to the energy differences of Eq. (f20l) . Correlated and uncorrelated single 
exponential fits to the correlation functions are performed to determine the energy of the 
two B-mesons, and the Jackknife method is used to determine the uncertainty. Further, 
Jackknife is used to determine the correlated difference in energy between the two .B-mesons, 
and twice the single 5-meson mass. 
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r 


A(r, 16) 


^(i-,16,6- 1 ) 


A(r, oo) 


B{r,oo,b- 1 ) 


(1,0,0) 
(1,1,0) 
(2,0,0) 
(2,1,0) 
(3,0,0) 
(A 0^1 

(5,0,0) 
(6,0,0) 
(7,0,0) 
(8,0,0) 


+0.2656 
+0.2011 
+0.1203 
+0.1144 
+0.1206 
—0 01 7fi 
+0.0592 
-0.1055 
+0.0409 
-0.1593 


-0.04261 
-0.1791 
-0.2352 
-0.2130 
-0.0750 
—0 1Q75 
+0.0874 
-0.1359 
+0.2310 
-0.1247 


0.2654 
0.2102 
0.1629 
0.1370 
0.1083 

0.0370 
0.0176 
0.0054 
0.0020 


-0.0938 
-0.2236 
-0.2451 
-0.2750 
-0.2357 
-0 9048 

-0.1653 
-0.1261 
-0.0921 
-0.0647 


TABLE V: The functions A(r, L) and B(r , L, ji) that contribute to the difference between the 
finite-lattice spacing and continuum OGE potentials at finite-volume. 


r 


6g*(f,16) 




3;NLO K > 




(1,0,0) 
(1,1,0) 
(2,0,0) 
(2,1,0) 
(3,0,0) 
(4,0,0) 
(5,0,0) 
(6,0,0) 
(7,0,0) 
(8,0,0) 


0.9229 
0.6406 
0.3762 
0.3941 
0.7327 
269.6 
2.092 
1.904 
16.87 
1.479 


-0.0794(40) 
-0.0425(21) 
-0.01798(90) 
-0.01529(77) 
-0.01202(60) 
+0.001319(66) 
-0.00354(18) 
+0.00525(26) 
-0.001747(87) 
+0.00595(30) 


-0.0828(41) 
-0.0524(26) 
-0.0272(14) 
-0.0227(11) 
-0.01398(70) 
-0.00254(13) 
-0.00217(10) 
+0.00348(17) 
+0.000836(42) 
+0.00473(24) 


-163.8(8.6) 
-103.7(5.4) 
-53.8(2.8) 
-45.0(2.4) 
-27.7(1.4) 
-5.04(26) 
-4.30(23) 
+6.90(36) 
+1.654(87) 
+9.37(49) 



TABLE VI: Corrections to the potential between two static color 3 sources combined into the 
3 representation, displaced by f computed on 16 3 x 32 DBW2 lattices with a lattice spacing of 
b = 0.0997 ± 0.0015 fm. The quoted uncertainty is due to the uncertainty in a s (2 GeV), and that 
of the lattice spacing. The NLO result does not use the BLM scale, q*, but the sum of the LO and 
NLO contributions. 



The effective mass plots for the correlators defining the central potentials at |r| = 
0,1,2,3,4,5,6,7,8 are shown in Fig. H] and those for the potentials at displacements of 
f = (1, 1, 0) and (2, 1, 0) are shown in Fig. [5l It is not possible to further decompose these 
latter potentials into the central and tensor components, without additional information. 
However, given that the tensor potentials at all the other displacements are found to be 
very small, it is not unreasonable to assume that they are also small for these displace- 
ments, and therefore we can assume that they provide a good determination of the central 
potentials alone. For a number of combinations of /, si and f , it was not possible to extract 
a signal and these points are omitted in the effective mass plots and tables below. 
After applying the perturbative one-loop matching discussed above, we determine the vari- 
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FIG. 4: The effective mass plots for the central potentials at |r| = 0,1,2,3,4,5,6,7,8 for each 
spin-isospin channel. Red stars correspond to = (0,0), green squares to (I,si) = (1,0), 

blue triangles to (/, si) = (0,1) and magenta diamonds to (I,si) = (1,1). Extracted masses and 
uncertainties are shown as the shaded regions in channels where a signal can be extracted, and are 
given in Tables rvTTllXl 



ous finite-volume potentials. The central potentials extracted from the lattice calculation in 
each of the spin-isospin channels are given in Tables IVIIllXl The lattice central potentials, 
V/>f > are shown in Fig. [6], and the central potentials with the leading order finite-lattice 

spacing correction included (as discussed above), Vjg , are shown in Fig. [71 In each channel 
there is a clean signal for the central potentials, with two or more of the displacements 
having potentials that are clearly non-zero. For the si = 1 channels the tensor potentials 
are found to be consistent with zero and are smaller than ~ 40 MeV in both cases for 
the entire range of displacements (the tensor potentials were also found to small and poorly 
determined in Ref. 16|). At large distances, this is not consistent with our expectations from 
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FIG. 5: The effective mass plots for the potentials at displacements f = (1, 1, 0) and (2, 1, 0) which 
are linear combinations of the central and tensor potentials. Red stars correspond to (I, s/) = 
(0,0), green squares to (I, si) = (1,0), blue triangles to (I, si) = (0,1) and magenta diamonds to 
(I, si) = (1,1). Extracted masses and uncertainties are shown as the shaded regions in channels 
where a signal can be extracted, and are given in Tables IVIIIIXl 



f 


b£§> B (r) 






Vjtl =0 (r) [MeV] 





0.783(11)(02) 


-0.325(13) (10) 


— oo 


— oo 


1 


0.960(12)(08) 


-0.148(14)(13) 


-0.231(15)(13) 


-456(30)(25) 


V2 


0.990(27)(07) 


-0.118(28)(12) 


-0.170(28)(12) 


-337(56)(24) 


2 

V5 
3 


1.056(13)(08) 


-0.052(15)(13) 


-0.079(15)(13) 


-156(30)(25) 


1.097(20)(15) 


-0.011(21)(18) 


-0.025(21) (18) 


-49(42)(36) 


4 


1.132(12)(06) 


0.024(14)(12) 


0.022(14)(11) 


43(28)(23) 


5 


1.143(09)(05) 


0.035(12)(11) 


0.033(12)(11) 


65(23)(22) 


6 
7 


1.150(09)(05) 


0.042(12)(11) 


0.046(12)(11) 


90(23)(22) 


8 











TABLE VII: The extracted lattice and continuum potentials in the I = S = channel. The first 
uncertainty is statistical, while the second is systematic. 

the NN system at the physical value of m w , but may result from the relevant B^'B^M 
couplings (M represents the various mesons) being small or giving rise to cancellations, or 
from the unphysically large pion mass. Further studies of this issue are warranted. The 
measurements at f = contain no information about the continuum potentials, which di- 
verge as ~ a(r _1 )/r. As discussed in Sections IIHI and HVl the lattice energies measured 
for coincident 5-mesons, in fact, determine the energies of the A& and E& and their exotic 
partners. 
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f 


6£g*(f) 






Vj% s=0 (f) [MeV] 




1 

V2 
2 

V5 
3 

A 


1.278(10)(05) 
1.199(07)(03) 
1.184(08)(01) 
1.158(07)(02) 
1.158(09)(01) 
1.142(06)(06) 


+0.170(12)(11) 
+0.091(10)(10) 
+0.076(11)(10) 
+0.050(10)(10) 
+0.050(12)(10) 
+0.034(10)(12) 


+oo 
+0.133(10)(10) 
+0.102(11)(10) 
+0.064(10)(10) 
+0.062(12)(10) 
+0.041(10)(12) 


+oo 
+262(21)(21) 
+203(22)(20) 
+126(20)(20) 
+122(23)(20) 
+82(19)(23) 


5 
6 
7 


1.134(05)(04) 
1.140(07)(04) 


+0.026(09)(11) 
+0.032(10)(11) 


+0.027(09)(11) 
+0.030(10)(11) 


+54(18)(21) 
+60(20)(21) 


8 


1.138(07)(03) 


+0.030(10)(10) 


+0.028(10)(10) 


-\-Ou[ZU)[Zl ) 


TABLE VIII: The extracted lattice and continuum potentials in the I = 
first uncertainty is statistical, while the second is systematic. 


1,5 = channel. The 


f 


b£* B (r) 


bvjT s ( =l(r) 


bvff s=1 (r) 


m 

Vjr l S =1 (f) [MeV] 



1 


0.850(06)(04) 
1.000(12)(18) 


-0.258(09)(11) 
-0.108(14)(21) 


— oo 
-0.191(15)(21) 


— oo 
-377(30) (41) 


2 










3 
4 
5 
6 
7 


1.128(18)(08) 
1.140(10)(10) 
1.150(10)(10) 
1.145(10)(15) 
1.1544(12)(09) 
1.149(10)(09) 


+0.020(20)(13) 
+0.032(12)(14) 
+0.042(12)(14) 
+0.037(12)(18) 
+0.036(14)(13) 
0.041(12)(13) 


-0.003(20)(13) 
+0.018(12)(14) 
+0.040(12)(14) 
+0.035(12)(18) 
+0.040(14)(13) 
0.042(12)(13) 


-05(38)(25) 
-\-^(\(9 K ,)(9R) 

+78(25)(28) 
+69(25)(36) 
+79(28) (27) 
+83(25)(27) 


8 











TABLE IX: The extracted lattice and continuum central potentials in the 7 = 5 = 1 channel. The 
first uncertainty is statistical, while the second is systematic. The asterisks attached to the r = y/2 
and y/5 potentials indicate that these are the sum of the central and tensor contributions. 

C. Potentials with t-Channel Quantum Numbers 

Up until this point we have classified the potentials between the B-mesons in terms of 
the s-channel quantum numbers, the total isospin and spin of the ldof. These potentials are 
extracted from the energies calculated on the lattice by subtracting twice the £>-meson mass. 
The statistical and systematic uncertainty in determining the B-meson mass propagates 
through to all four central potentials, leading to larger uncertainties in the potential than 
from the calculation of the energy of the two .B-mesons alone. Motivated by the success of 
traditional nuclear physics phenomenological potentials constructed from the exchange of 
mesons in the t-channel, we have formed linear combinations of the s-channel potentials to 
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r 








vgJ^Cf) [Me 





1.253(08)(02) 


+0.145(11)(10) 


+oo 


+oo 


1 


1.172(07)(05) 


+0.064(10)(11) 


+0.106(10)(11) 


+209(21)(22) 


V2* 
2 


1.163(06)(01) 


+0.055(10)(10) 


0.081(10)(10) 


+161(19)(20) 


V5* 
3 


— 

1.136(08)(03) 


— 

+0.028(11) (10) 


— 

+0.040(11) (10) 


— 

+78(22)(21) 


4 


— 

1.128(07)(04) 


— 

+0.020(10) (11) 


— 

+0.021(10)(11) 


— 

+42(20)(21) 


5 


1.123(07)(06) 


+0.015(10)(11) 


+0.016(10)(12) 


+32(20)(23) 


6 


1.134(06)(02) 


+0.026(10)(10) 


+0.024(10) (10) 


+48(19)(20) 


7 


1.149(10)(09) 


0.041(12)(13) 


0.042(12)(13) 


+83(25)(27) 


8 


1.131(06)(03) 


0.023(10)(10) 


0.021(10)(10) 


+41(19)(21) 



TABLE X: The extracted lattice and continuum central potentials in the 1 = 0,5' = 1 channel. 
The first uncertainty is statistical, while the second is systematic. The asterisks attached to the 
r = y/2 and v5 potentials indicate that these are the sum of the central and tensor contributions. 

give potentials with well-defined spin and isospin quantum numbers that can be identified 
with the exchange of one or more hadrons. The central potential can be decomposed as 

V I , Sl (\r\) = V 1 (\r\) + a x -a 2 V a {\r\) + a x ■ a 2 n ■ r 2 V ffT (|r|) + n • r 2 V T (\r\) , (30) 

where it is straightforward to show that, in terms of the s-channel central potentials, 

Vi = 4 (^o,o + 3Vo,i + 3^,0 + 9^,0 , 
lb 

V = — {-V ofi + - 3T4 )0 + 3V ltl ) , 
lo 

lb 

V T = ^ {-V 0l0 - 3V ,i + Ki,o + 3K,i) • (31) 

An important point to observe is that the three potentials, V a , V UT and V T can be extracted 
from the lattice calculation without reference to the 5-meson mass. This is not true for V\. 
It was pointed out in Ref. [24| that the rj' double-pole that is present in quenched calculations 
will dominate the interactions between nucleons and I?-mesons at long-distances. However, 
we see that this can only contribute to V a , and not to the other three potentials. Further, 
the exchange of a single 7r will contribute only to the V aT potential. Therefore, the potential 
V T does not receive contributions from hairpins nor from 7r-exchange, and does not depend 
upon the 5-meson mass extraction from the lattice calculation. It is expected to be clean, 
and determined by short-range and medium range interactions. 

The finite- volume potentials calculated on the lattice, V^^i(f), are given in Tables IXll- 
IXIVl and are shown in Fig. [HJ Effective mass plots for these potentials are shown in Figs. 
M and (TDJ The potentials corrected for the finite-lattice spacing contributions to OGE be- 
tween the heavy quarks are also given in Tables IXItiXIVl and are shown in Fig. [TTJ The 
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FIG. 6: The finite- volume and finite-lattice spacing central potentials, , extracted from 

the lattice calculation. The shaded regions are simple fits to guide the eye. Red stars, green 
squares, blue triangles and magenta diamonds correspond to (I,si) = (0,0), (1,0), (0,1), (1,1) 
respectively. 



f 


bV^ tm (r) 


bV^ L \f) 


v£ L \r) [MeV] 





-0.050(17) (06) 


— oo 


— oo 


1 


-0.0174(33) (08) 


-0.0337(34) (08) 


-66.8(6.8)(1.7) 


V2* 
2 


-0.0081 (25) (02) 


-0.0178(25) (02) 


-35.1(5.1)(0.3) 




-0.0016(25) (09) 


-0.0049(25)(09) 


-9.7(5.0)(1.8) 


3 


+0.0035(13) (00) 


+0.0008(12)(00) 


+1.6(2.6)(0.1) 


4 


+0.0048(13) (08) 


+0.0035(13)(08) 


+6.9(2.6)(1.5) 


5 
6 


+0.0054(12)(05) 


+0.0044(12)(05) 


+8.8(2.4)(1.0) 


7 


+0.0054(18)(02) 


+0.0054(18)(02) 


+10.7(3.6)(0.3) 


8 


+0.0047(16)(03) 


+0.0053(16)(03) 


+10.5(3.1)(0.6) 



TABLE XI: The extracted lattice and continuum central potential , defined in Eq. (|31l) . The 
first uncertainty is statistical, while the second is systematic. The asterisks attached to the r = v2 
and y/b potentials indicate that these are the sum of the central and tensor contributions. 



uncertainties in these potentials are seen to be significantly smaller than those of Vi iSv in 
part due to the fact that the 5-meson mass extraction, and its associated uncertainty, does 
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FIG. 7: The finite-volume central potentials extracted from the lattice calculation including 
the leading lattice correction to OGE, Vj^. The shaded regions are simple fits to guide the 
eye. Red stars, green squares, blue triangles and magenta diamonds correspond to (I, si) = 
(0,0), (1,0), (0,1), (1,1) respectively. 



r 




bv£\r) 


Vjr\r) [MeV] 





-0.0539(56)(07) 


-co 


— oo 


1 


-0.0228(24)(08) 


-0.0383(25) (08) 


-75.9(5.2)(1.5) 


V2* 


-0.0158(21)(02) 


-0.0256(22)(02) 


-50.8(4.4) (0.3) 


2 


-0.0064(22)(07) 


-0.0115(22) (07) 


-22.7(4.3)(1.3) 


V5* 


-0.0062(42)(05) 


-0.0105(42) (05) 


-20.7(8.4) (1.0) 


3 


-0.0003(23) (12) 


-0.0029(22) (12) 


-5.7(4.5) (2.3) 


4 


+0.00205(82) (50) 


+0.00158(82)(50) 


+3.1(1.6)(1.0) 


5 
6 


+0.00343(09) (32) 


+0.0030(09)(03) 


+6.0(1.9)(0.6) 


7 

8 


+0.00294(55)(07) 


+0.00309(55)(07) 


+6.1(1. 1)(0.1) 



TABLE XII: The extracted lattice and continuum central potential V& T , defined in Eq. (|3ip . The 
first uncertainty is statistical, while the second is systematic. The asterisks attached to the r = v2 
and y/b potentials indicate that these are the sum of the central and tensor contributions. 

not contribute. The most striking potential is V T ; it is clear that this potential is of shorter 
range than V a and V aT , due to the absence of OPE and OHE. 
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f 


b V T (r) 


o V-r (r) 


Vt \r) [MeVJ 





-0.0390(49) (03) 


— oo 


— oo 


1 


-0.0136(39)(02) 


-0.0292(40)(02) 


-57.7(7.8)(0.4) 


/a* 

V2 


-0.0033(31)(05) 


-0.0131(31)(05) 


-25.9(6.3)(l.O) 


2 


+0.00539(66)(02) 


+0.00029(70) (02) 


+0.58(l.4)(0.0) 




+0.00587(53)(15) 


+0.00161(57)(15) 


+3.2(l.2)(0.3) 


3 


+0.00787(71) (12) 


+0.00525(73)(l2) 


+10.4(1.4) (0.2) 


4 


+0.00767(61)(22) 


+0.00720(6l)(22) 


+14.2(1.2) (0.4) 


5 
6 


+0.00664(97) (04) 


+0.00623(97) (04) 


+12.3(1.9)(0.1) 


7 


+0.00509(50) (07) 


+0.00525(50)(07) 


+10.4(1.0)(0.1) 


8 


+0.00455(62) (24) 


+0.00544(62) (24) 


+10.8(1.2)(0.5) 



TABLE XIII: The extracted lattice and continuum central potential Vt , denned in Eq. (I3ip . The 
first uncertainty is statistical, while the second is systematic. The asterisks attached to the r = v2 
and y/5 potentials indicate that these are the sum of the central and tensor contributions. 



r 


bVi att{L \r) 


b v} L) (r) 


v} L \r) [MeV] 





-0.117(16) (03) 


— oo 


— oo 


1 


-0.0616(96) (25) 


-0.0520(96)(25) 


-103(19)(04) 


V2* 


-0.0193(74) (04) 


-0.0210(74) (04) 


-42(15)(0.7) 


2 


-0.0135(58) (18) 


-0.0101(58)(18) 


-20(12) (4) 




-0.0020(78) (10) 


-0.0019(78) (10) 


-4(16)(2) 


3 


+0.0030(43)(17) 


+0.0037(43)(17) 


+7.4(8.4) (3.3) 


4 


+0.0093(45)(10) 


+0.0094(45)(10) 


+16.6(8.9)(2.0) 


5 
6 


+0.0069(19) (01) 


+0.0064(19) (01) 


+12.6(3.8)(0.1) 


7 


+0.0074(43)(14) 


+0.0057(43)(14) 


+11.2(8.4)(2.8) 


8 


+0.0051(29)(07) 


+0.0055(29)(07) 


+10.8(5.7)(1.4) 



TABLE XIV: The extracted lattice and continuum central potential , defined in Eq. (I3ip . The 
first uncertainty is statistical, while the second is systematic. The asterisks attached to the r = y/2 
and y/5 potentials indicate that these are the sum of the central and tensor contributions. 



D. Extrapolation to Infinite- Volume Potentials 

The final stage of analysis is to use the extracted finite-volume potentials in either the 
s- or t-channels to determine the infinite-volume forms. At short distances, |r| < A" 1 
(|f | ^2 for our analysis), the infinite- volume extrapolation must be done empirically, fitting 
functions with the correct long-distance behavior to the results of lattice calculations in 
multiple volumes. In principle, this extrapolation can be performed systematically for larger 
|r| as effective field theory describes the potential in this regime. Here we explore how the 
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FIG. 8: The central finite- volume lattice potentials V^^^\, as defined in Eq. (|3T1) . The statistical 
and systematic errors have been added in quadrature. 



matching to EFT can be implemented, focusing on the isovector potentials. 

In QCD, the long range pieces of the infinite-volume, t-channel isovector potentials are 
expected to have the form 



g 2 ml 

247T/2 



-m w |r| 



V^\\r\) + 



(32) 



V (oo), 



|r|- 



47T 



K (27r) (|r|) + 



(33) 
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where I/.^ 271 ") and V^? 7r ) are the two pion exchange potentials defined in Ref. [43j (with nucleon 
couplings replaced by the relevant couplings of the B sector), f n ~ 132 MeV, g is the chiral 
coupling of pions to heavy mesons occurring in the heavy meson chiral perturbation theory 

and g p is a phenomenological BBp coupling. The ellipses 



Lagrangian [4J, |45|, |46|, L 
denote contributions suppressed at large separations. 

Ideally, lattice determinations of the meson masses and potentials at long distances could 
be used to fit the couplings in the above equations (the two pion contributions contain ad- 
ditional parameters). However, a number of issues complicate this analysis. The quenched 
nature of our calculations introduces artefacts as in this case, the rj' meson remains degen- 



erate with the pions but has a modified propagator [H, 13 



G^q 2 



+ 



i(M 2 - a^q 2 



q 2 — m 2 + ie (q 2 — m 2 + ie) 2 



(34) 
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FIG. 9: The effective mass plots for the central potentials at |f| = 0,1,2,3,4,5,6,7,8 for each 
t-channel potential. Red stars correspond to V a , green squares to V aT , blue triangles to V T and 
magenta diamonds to V\. Extracted masses and uncertainties are shown as the shaded regions in 
channels where a signal can be extracted. 



(M and a$ are couplings occurring in the quenched chiral Lagrangian [49|, |5C|) that pro- 
duces unphysical components of the potential. In particular, both of the isovector t-channel 
potentials receive contributions from one-pion-one-7/ exchange that are longer range than 
the two pion exchange contributions. These contributions are calculable, but involve ad- 
ditional low energy constants. Additional issues are introduced by the unphysically large 
quark mass used in our calculations with the identification of the dominant contribution in 
V T depending on the quark mass. At the physical mass, single p exchange is sub-dominant 
to two pion exchange, however here, m p < 2m n so it is p-exchange that persists to the 
longest distance. In our calculation, 2m 7r ~ A x and, in both channels, the two Goldstone 
boson exchange contributions are indistinguishable from short distance contributions not 
describable in EFT. Finally, the formula for the potential at finite-volume in Eq. (1261) is 



valid only for single particle exchange and is significantly modified if two or more particle 
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FIG. 10: The effective mass plots for the potentials at displacements f = (1, 1,0) and (2, 1,0) Red 
stars correspond to V a , green squares to V aT , blue triangles to V T and magenta diamonds to V\. 
Extracted masses and uncertainties are shown as the shaded regions in channels where a signal can 
be extracted. 
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FIG. 11: The central finite-volume potentials V CjCTT)Tj i, as defined in Eq. (f3Tj) (the lattice potentials 
plus the leading lattice spacing corrections to OGE). The statistical and systematic errors have 
been added in quadrature. 



exchange effects are included at infinite- volume; the two particles can interact with sources 
in different periodic copies. 

Since only the longest range contribution to the potential in each channel can be identified, 
we fit our results at large separations, |r| > A" 1 , using the finite- volume versions (computed 
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FIG. 12: Fits to the finite- volume isovector i-channel potentials. The dashed lines correspond to 
the finite- volume fits to the lattice data, and the solid curves are the infinite- volume extrapolations. 



using Eq. 



of the simplified infinite-volume potentials, 



-nv|r| 
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+ a' 



e -A x |r| 
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-m p \r\ 



-Axl'l 



+ ^ 



(36) 



Using the measured values and uncertainties of and m p and the physical value of f n we 
first determine the couplings g and g p by setting a x = a! = and fitting the finite- volume 
potentials at the two largest separations. 9 These fits are shown by the dashed red curves in 
Fig. H2] and the resulting couplings are found to be 



g p = 2.17 ±0.08, 



g = 0.57 ±0.06. 



(37) 



These couplings are stable under decreasing the minimum separation toward the point where 
the finite- volume potential crosses zero however the \ 2 of the fit worsens. Having determined 
these parameters, we reconstruct the infinite-volume potentials that are shown in the figure 
as the solid red lines. 

If the couplings a W in Eqs. (1331) and fl3"B"l) are included as fit parameters, we obtain instead 



9 P 



3.02 ±0.09 



g = 0.69 ±0.03 



(38) 



with the finite-volume fits and their infinite-volume reconstructions shown as the dashed- 
and solid- blue curves in Fig. [12j In this case we have set the minimum separation, r m i n , 
used in our fits to be 2b for V aT and 36 for V T although the fits vary only slightly under 
changes of r m ; n from b to 56. Averaging the two sets of extractions, we find 



<7 P = 2.6 ±0.1 ±0.4, 



g = 0.63 ±0.05 ±0.06 



(39) 



Simple fits using the infinite- volume long range behaviour were considered in Ref. 16 1 
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where the second error is an estimate of systematic errors determined by differences between 
the two fits and variation of the fit range. These numbers represent our best estimates of 
the couplings but we caution that we are currently unable to investigate the full systematics 
of this determination. Further refinement would require lattice calculations at a range of 
different volumes, lattice spacings and quark masses. 

In both isovector channels, the agreement of the two infinite-volume extractions at large 
separations suggests that the long range piece of the extraction is robust. The pion cou- 
pling, g, is related to the forward limit of (B\ \B*), the matrix element of the isovector 
axial- vector current, through PC AC and the value we extract is consistent with direct de- 
terminations of the quenched axial coupling: 0.42(4)(8) Q, 0.69(18) 0, 0.48(3)(11) (£}, 
0.517(16) [Hi]]. We note that extraction of this coupling from the potential does not require 
renormalisation of axial current and suffers from different systematic effects. Agreement 
between the two procedures is encouraging. 

The isoscalar channels suffer from more severe unphysical artefacts in quenched QCD 
and we are not able to extract meaningful information from the long distance potentials. 
For V a , EFT predicts a long range single Goldstone boson exchange potential |24j 



a " ] ^\ 2^ 6 



(40) 



(go is the rf axial-coupling occurring in the quenched heavy- meson chiral Lagrangian [471.1481]) 
with a long distance exponential tail dominating. Unlike the other channels, the suppression 
of the sub-leading contribution is not exponential and our data is insufficient to resolve these 
pieces. V\ is not determined by single particle exchange (though many phenomenological 
approaches in the nucleon sector include exchange of the <r(550) resonance) and in this chan- 
nel our data are particularly poor. In both isoscalar channels two-?/ exchange is also present 
and enhanced compared to two-pion, and one-r/-one-pion- exchange, further polluting the 
signals. Larger volumes and multiple quark masses will be needed to perform extractions of 
couplings in these channels. 



VI. DISCUSSION 

We have studied the potentials between two 5-mesons in the heavy quark limit. The 
calculations were performed on 16 3 x 32 quenched lattices with a spatial length of ~ 1.6 fm, 
and with a quark mass such that ~ 400 MeV. The leading lattice space corrections to the 
one-gluon-exchange potential between the two heavy quark propagators in the finite- volume 
were included in order to extract the physical potential between I?-mesons in the continuum 
but at finite- volume. We find clear evidence of repulsion between the 5-mesons in the I ^ S; 
channels and attraction in the I — s\ channels. Three of the four potentials defined with 
i-channel spin-isospin quantum numbers have significantly smaller uncertainties than the 
potentials defined with s-channel quantum numbers. From the large separation behaviour 
of these potentials at finite- volume, S-meson couplings to the 7r and p were extracted. 

This calculation can be improved in a number of areas but shows that a rigorous first 
principles calculation of the S-meson potential is achievable in the near future. The next 
stage of our study will progress from unphysical quenched QCD to fully dynamical QCD. 
This is mandatory for connection to the real world but will also significantly simplify the 
analysis of the long range potential using EFT. To separate the different components of 
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the potential in the short-, intermediate- and long- range regimes requires multiple volumes 
and quark masses. Finally calculations at a number of different lattice spacing are required 
to control the remaining discretisation effects. Completion of this ambitious program will 
provide deep insight into the BB system and, ultimately, nuclei. 
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APPENDIX A: POTENTIALS FROM LATTICE WAVEFUNCTIONS 

Recently, it has been claimed that the nucleon-nucleon potential can be extracted from 
the lattice wavefunctions of two nucleons |56j, extending the technique that CP-PACS has 
successfully used to determine I = 2 tttt scattering parameters [5?]]. In this appendix, we 
question the validity of that calculation. 

An interpolating field for the nucleon is, 

dl(x,t)i = e abc & (q a > T C l5 r 2 q b )(x,t) , (Al) 

where i is a Dirac-index, a is an isospin index, and a, b, c are color indices. This operator 
has a non-zero overlap with a nucleon momentum-eigenstate 

(01^(0, 0)^(p)) = ZsWSapS* . (A2) 

The two-nucleon correlation function measured on the lattice in Ref. 56[ is 

GW(x,y,t) = (oio^tyMy^ypmo} 



g—E„t 



= £ (0|d 1 (x,0)|,d 1 (y,0)j|V„>^n|J(0)|0> — , (A3) 

where J is a wall-source on the initial time-slice to = 0, and \ip n ) are the eigenstates of the 
Hamiltonian in the finite- volume. In particular, \ip n ) are states of definite baryon number 
(B = 2), isospin and transformation under the hyper-cubic group. Setting (ip n \J(to)\0) = 
A n (to), at long times the correlation function becomes 

P -E t 

GW(x,y,t) - 4,(0) (0iO 1 (x,0) l a O 1 (y, 0)^ > -r=- , (A4) 
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and Ref. |56| asserts that (0|Oi(x, t^O^y, t ) J p\A) = (0|iV(x, t Y a N(y, t )^|Vo> = ^ 
is the Bethe-Salpeter wavefunction (their Eq. (4)). From this definition, they generate a 
nucleon-nucleon potential via 

V{r) = E + = E + , (A5) 

where $ at a? is the projection of the Bethe-Salpeter wavefunction onto definite isospin and 
transformation under the hyper-cubic group, fi is the reduced mass of the two-nucleon 
system. However, this identification of the Bethe-Salpeter wavefunction is incorrect and the 
most general form for the matrix element is 

(0|a i (x,t )id 1 (y,to)j|^> = 4 S ^(ir|) (0\N(x,t o y a N(y,t %\A) + ■ ■ ■ , (A6) 
( s 1} 

where Z N ' N (|r|) is an unknown function that depends upon the spin, isospin and structure 
of the composite sink, Oi(x, t ) l a Oi(y, t )p, and where r = x — y. The ellipsis denotes 
additional contributions from the tower of states of the same global quantum numbers. 
With this complete form of the matrix element, it is not possible to determine the potential 
from Gnn without additional information. In the limit |r| — ► oo, Z^n (| r |) ~~ > 1 and the 
additional terms in Eq. (IA6I) containing p > 2 particles are suppressed. Consequently the 
scattering parameters can be rigorously extracted as has been done for the case of the 



pion [57 . 



APPENDIX B: FINITE LATTICE SPACING CORRECTION TO THE POTEN- 
TIAL 

The finite lattice spacing correction to the potential (in the 3 color channel) is given by, 



SV (L) -(r) 



~37T 2 6 



2^ 3 

T 



' L\ e *27rn-f/L 

2^1 ? In! 2 



lml<£ 



E el2 ™' f/L G 00 (n x ,n y ,n z ,0) 



(Bl) 



arising from the difference between continuum and lattice one gluon exchange evaluated at 
finite- volume. 

The lattice contribution to this expression is simple to evaluate for the DBW2 action 
(the full form of the improved gluon propagator is given in Ref. 35j), however calculating 
the continuum contribution is somewhat subtle. Difficulties arise in both the infrared and 
ultraviolet regimes. Both the lattice and continuum finite- volume sums are IR divergent, 
however provided both are regulated in the same way a sensible result ensues; the simplest 
procedure is to omit the zero-mode 10 . 



10 Any regularization is equally valid as the mode expansion of the perturbative gluon propagator is in- 
trinsically ill-defined in the IR region. Differences in IR regularization lead to 0(b 2 ) differences in the 
perturbative corrections to the potentials, parametrically smaller than the effects of the Wilson fermion 
discretisation used herein. 
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While the continuum contributions to A and B, defined in Eq. (1291) . are strictly UV con- 
vergent, that convergence is highly oscillatory. Computing these contributions is simplified 
by the use of the Poisson summation formula which allows the sum to be rewritten as 

„27rm-f/L p 2nin-r/L „27rin-f/L 

\ — 2 \ — T \ — T 

^o^T = X n Vo|n| 2 (|n| 2 + ^) + n ^ (|nP + ^) 

„27rin-f/L i r „— 2n\m+r/L\x 

2 V - - + — e ~ 2nx W /L + 7T V - (B2) 



which is independent of the value of x. The sums on the rhs of this expression are more 
convergent than that on the lhs and can be numerically evaluated reliably. Similar techniques 
allowed us to deal with the analogous differences defining the function B. 

In the limit that |f| — >• 0, the continuum contribution is singular, leading to a correction 
factor of 

<» - + - • (B3) 

which is nothing other than the strong Coulomb interaction between the heavy quarks in 
the continuum. In the continuum limit and infinite-volume limit, b <C r <C L, the leading 
correction factor is found to be 

*Or)- -5^(1 + 12* -124)+..., (B4) 

where we have used the BLM procedure to set the scale. This improved perturbative shift 
can be eliminated for suitable choices of C\. Clearly, the Luscher-Weisz-improved value of 
C\ = Y2 maximally improves the lattice calculation. That is to say that, neglecting the small 
cf contribution, the correction factor that must be applied to the lattice calculation in order 
to recover the continuum potentials is minimized by Luscher-Weisz-improvement. 
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